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Abstract 

We report a detailed examination the magnetic structure of anhydrous cupric chloride CuCl2 
carried out by powder neutron diffraction, magnetic susceptibility and specific heat measurements 
on polycrystalline and single crystal samples as well as an evaluation of the spin exchange inter- 
actions by first principles density functional theory (DFT) calculations. Anhydrous CuCl2 shows 
one dimensional antiferromagnetic behavior and long range antiferromagnetic ordering below a 
Neel temperature of 23.9 K. Neutron powder and single crystal diffraction reveal that, below 23.9 
K, CuCl2 undergoes a phase transition into an incommensurate magnetic structure (propagation 
vector (1,0.2257,0.5) with a spin-spiral propagating along h and the moments confined in the he 
crystallographic plane. Our DFT calculations show that the spin-spiral results from competing 
ferromagnetic nearest neighbor and antiferromagnetic next-nearest neighbor spin-exchange inter- 
action along the spin chains. Implications for possible multiferroic behavior of CuCl2 are discussed. 
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I. INTRODUCTION 



Cupric chloride CuCl2 has a crystal structure in which layers of composition CuCl2 stack 
along the c-direction with van der Waals interactions between them. Each CuCl2 layer con- 
sists of CuCl2 ribbons made up of edge-sharing CUCI4 square planes such that the axial 
positions of each CuCU square are capped by the CI atoms of its adjacent ribbons with long 
Cu...Cl distances. Thus, each Cu^"^ (d^, S = 1/2) ion is contained in an axially elongated 
CuCle octahedron, as expected for a Jahn- Teller active Cu^"*" ion. The magnetic suscep- 
tibility measurements exhibit a broad maximum, characteristic of a short-range ordering 
predicted for low-dimensional magnetic systems.->2i>^'^ Historically, CUCI2 was among the 
very first magnetic systems whose bulk magnetic properties were analyzed in terms of a lin- 
ear spin chain model. Heat capacity and NQR measurements indicated a Neel temperature 
Tn ~ 23.9 K for CUCI2. To the best of our knowledge the magnetic structure of CUCI2 has 
remained unresolved until now. 

CuCl2 is a chemically simple one-dimensional (ID) AFM spin 1/2 quantum chain sys- 
tem. However, by analogy to the CUO2 ribbon chains found in the magnetic oxides such as 
LiCuV04 and LiCu202, the next-nearest-neighbor (NNN) spin exchange coupling along the 
chain may be essential in addition to the nearest-neighbor (NN) spin exchange in under- 
standing the magnetic properties of CuCl2. For these CUO2 chains, made up of edge-sharing 
CUO4 square planes, the NN spin exchange is ferromagnetic (FM), while the NNN spin 
exchange is strongly AFM, and the resulting geometric spin frustration gives rise to an in- 
commensurate spiral spin ordering and ferroelectric polarization. Similarly, CUCI2 might 
undergo a spiral spin ordering within each Cu chain and a long-range Neel ordering as do 
LiCuV04 and LiCu202. In the present work we examine the magnetic structure of CuCl2 in 
detail by powder neutron diffraction, magnetic susceptibility and specific heat measurements 
as well as by evaluating its spin exchange interactions on the basis of first principles density 
functional DFT calculations. 

II. EXPERIMENTAL 

Single crystals of CuCl2 were grown in quartz glass ampoules by the Bridgman technique 
using commercially available powder (Alfa Aesar, ultra dry, purity 99.995%). Special care 
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was taken in the design of the ampoules to support the elevated vapor pressure of chlorine 
above the melting point of CuCl2. Well-crystallized needle-shape crystals, obtained with 
the needle axis oriented along the crystallographic b axis were found to be twinned in the 
ac-plane. Due to the moisture sensitivity of CuCl2 all handling of the samples was done in 
an Ar filled glove box. Temperature and field dependent magnetizations (0 T < B^xt < 7 T) 
were measured using a superconducting quantum interference device magnetometer (MPMS, 
Quantum Design, 6325 Lusk Boulevard, San Diego, CA.). Magnetizations up to 50 T were 
measured in pulsed magnetic fields at the Laboratoire National des Champs Magnetiques 
Pulses (Toulouse, France) at temperatures of 1.4 and 4.2 K with the field aligned along 
the crystal b axis. The heat capacities of crystals (m ~ 20 mg) were measured using a 
commercial Physical Property Measurements System calorimeter (Quantum Design, 6325 
Lusk Boulevard, San Diego, CA.) employing the relaxation method in magnetic fields up to 
9 T aligned parallel and perpendicular to the crystal b axis. To thermally anchor the crystals 
to the heat capacity platform, a minute amount of Apiezon N grease was used. The heat 
capacity of the platform and the grease was determined in a separate run and subtracted from 
the total heat capacity. The heat capacity of powder samples (m ~ 300 mg) was measured 



in a home-built fu 
in detail in Refs. [s 



automatic Nernst-type adiabatic calorimeter similar to that described 
|. The powder samples were encapsulated in Duran glass flasks under ~ 
900 mbar ^He atmosphere to enable rapid thermal equilibration. The flasks were attached 
with a minute amount of Apiezon N vacuum grease to the sapphire platform which carries an 
deposited thin-film stainless steel heater and a calibrated Cernox CX-1050-SD thermometer 
(Lake Shore Cryotronics, Inc. 575 McCorkle Blvd. Westerville OH 43082). The heat 
capacities of the sapphire sample platform, the glass flask and the Apiezon vacuum grease 
were determined in separate runs and subtracted from the total heat capacity. EPR spectra 
were measured with a Bruker ER040XK microwave X-band spectrometer and a Bruker BE25 
magnet equipped with a BH15 field controller calibrated against DPPH. A stepper-motor- 
controlled goniometer was used to rotate the crystals around an axis perpendicular to the 
magnetic field plane. The crystals were oriented optically with a precision of ~ ±5°. 100 kHz 
field modulation was used to record the first derivative of the signal intensity. The spectra 
were fitted to a resonance curve with Lorentzian lineshape (zero dispersion) by varying the 
resonance position as well as the linewidth, intensity and background parameters. 

Powder neutron diffraction patterns were collected on the high-intensity two-axis diffrac- 
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tometer D20 (Institut Laue-Langevin, Grenoble) at a wavelength of 2.4 A and 1.88 A in 
the temperature range 1.8 K< T < 50K.- Single crystal diffraction was carried out on the 
four-circle diffractometer DIO (Institut Laue-Langevin, Grenoble) at a wavelength of 2.36 A 
on a crystal of ~50 mm^.- 

III. EXPERIMENTAL RESULTS 
A. Crystal Structure 

At room temperature, anhydrous cupric chloride CuCl2 was reported to crystallize with a 
monoclinic structure (space group C2 / m) with 2 formula units per unit celL-^^^ The crystal 
structure of CuCl2 consists of CuCl2 slabs parallel to the ac-plane that are made up of CuCl2 
chains running along the crystallographic h axis. Such slabs are interconnected along the c 
axis via van der Waals contacts (cf. Figure [1]). 

To confirm the crystal structure details of CuCl2 and to search for structural phase 
transitions at low temperatures, we collected a series of neutron diffraction patterns on 
the diffractometer D20 at temperatures 2 K < T < 50K. Patterns with sufficient intensity 
but improved resolution were collected using neutrons with a wavelength of A = 1.88 A 
(Ge monochromator) suitable for Rietveld profile refinements. These were done with the 
FuUProf program package^. The patterns were refined assuming the space group C2/m 
with the Cu atoms in the Wyckoff position (2a) (0,0,0) and the CI atoms in the position (4i) 
(x,0,2;) using a slightly asymmetric pseudo-Voigt profile shape to model the Bragg reflections. 
Within experimental error the atom coordinates x and z for the CI atoms are independent 
of the temperature and amount to x=0. 5065(2) and ^=0. 2365(3). Figure [2] displays the 
experimental pattern collected at 42.9 K and the results of the Rietveld refinement (A=1.88 
A ). Apart from changes in the lattice parameters, the structure remains unchanged with 
respect to the room temperature crystal structure of CuCl2 described by Wells and Burns 
et al.M^ The atom position parameters show no noticeable difference to those of the room 
temperature structure. The structure parameters at 2.5 and 42.9 K are summarized in Table 

m 

The temperature dependence of the lattice parameters is displayed in Figure [3l The 
lattice parameters show a monotonic decrease with temperature. Below ~ 15 K they level 
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FIG. 1: (Color online 



Crystal structure of CuCl2 projected along [010] drawn using parameters 



according to Refs. I9ljl0l. A unit cell is outlined by the black solid lines. The Cu atoms are 
represented by small (blue), and the CI atoms by large (orange) circles. The lower part of the 
figure highlights a CuCl2 ribbons with the h axis along the ribbons. The planes of the ribbons are 
parallel to the he plane. 

off and remain constant to the lowest temperatures. Near ~ 24 K, especially c and /?, exhibit 
noticeable changes in the slope, which indicate a magnetoelastic response of the lattice to the 
onset of long-range AFM ordering, as is evidenced by our magnetization and heat capacity 
experiments discussed in the following paragraphs. 

Single- crystal x-ray diffraction on crystals of CuCl2 were carried out using a laboratory 
based diffractometer (PDS II, Stoe & Cie., Darmstadt, Germany). These measurements 
confirmed the crystal structure determined from neutron powder diffraction and showed the 
crystals to be twinned in the ac- plane. The twinning matrix T relating the Miller indices 
of the twinned individua was determined to 
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FIG. 2: (Color online) Neutron diffraction pattern (red o) (A=1.888 A ) of a powder sample of 
CuCl2 collected at 42.9 K at the high intensity diffractometer D20 plotted in comparison with the 
pattern obtained from a profile refinement (solid black line). The difference is shown in the lower 
part by the (blue) solid line. The positions of the Bragg reflections used to simulate the patterns 
are indicated by the small (green) vertical bars in the lower part of the figure. 



fi -A 

r= -1 . (1) 

The twinning reverses the 6-axis and rotates the the axes in the reciprocal space, such 
that twin-related axes enclose an angle of 180° - /3 ~ 57°. 



IV. MAGNETIC PROPERTIES 



A. Magnetic Susceptibility Measurements 

Figure m displays the magnetic susceptibility of a polycrystalline sample of CuCl2. The 
susceptibility is characterized by a broad maximum centered at ~ 70 K and a Curie- Weiss 
law (Eq. behavior at high temperatures (see the detailed discussion of the susceptibility 
measurement on a crystal below) with a negative Curie- Weiss temperature of ~ -60 K and 
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T(K) a(A ) 6 (A ) c (A ) /? (°) method 

2.55 6.7986(3) 3.29418(12) 6.7718(3) 122.45(2) neut 

42.9 6.8010(3) 3.29437(12) 6.7768(3) 122.47(2) neut 

293 6.8973(4) 3.2961(2) 6.8160(4) 122.239(3) x-ray 

293 6.9038(9) 3.2995(4) 6.824(1) 122.197(8) x-ray (Ref. 10) 

TABLE I: Lattice parameters of CuCl2 obtained from neutron powder and x-ray diffraction patterns 
(Cu Kai radiation) at the indicated temperatures. Literature values are given for comparison. 

a powder gf-factor (see Eq. ([5]) of ~2.13. 

The negative Curie- Weiss temperature proves predominant AFM exchange interactions, 
and the broad maximum reveals AFM short range ordering indicative of a low-dimensional 
magnetic behavior. A kink-like change of the slope near ~ 24 K is due to the onset of long- 
range AFM ordering. A A-type anomaly is clearly revealed by the temperature derivative 



of the quantity Xmoi x T ('Fisher's heat capacity', cf. Ref. Il2l ) shown in the inset of Figure 
m from which the Neel temperature Tn is determined to be 23.9(3) K. coincides with 
the anomalies observed in the temperature-dependence of the lattice parameters (Figure [3]). 
When fitting a Curie- Weiss law to the high temperature data of the powder susceptibility 
it is found that the fitted parameters (the slope of the reciprocal susceptibility i.e. the 
effective magnetic moment, and the intercept) are strongly correlated and depend essentially 
on temperature-independent contributions to the susceptibility (diamagnetic and Van Vleck 
susceptibilities). 

In order to obtain a reliable fit of the magnetic susceptibility of CuCl2 over a wide 
temperature range and as well as the temperature independent diamagnetic and paramag- 
netic corrections to the susceptibility we determined the magnetic susceptibility on a crystal 
(5cxt=lT II b) up to 600 K. At high temperatures the magnetic susceptibility follows a 
Curie- Weiss law as demonstrated by the plot of the reciprocal susceptibility in Figure O 
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FIG. 3: Temperature dependence of the refined lattice parameters of CuCl2 resulting from Rietveld 
profile refinements of the neutron powder diffraction patterns (details see text). 
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FIG. 4: (Color online) Temperature dependence of the magnetic susceptibility (not corrected for 
temperature independent diamagnetic or Van Vleck paramagnetic contributions) of a polycrys- 
talline sample of CuCl2 determined in an external magnetic field of 1 T. The upper inset shows the 
reciprocal susceptibility with the (red) solid line indicating a Curie- Weiss law with an intercept on 
the abscissa at ^ -60K. The lower inset shows d/dT{xmoi x T ('Fishers heat capacity'). A sharp 
A-type anomaly at 23.9(3) K indicates the onset of long-range AFM ordering. 

Deviations from the Curie- Weiss law are seen already at temperature below ~ 175 K. 

At high temperatures, besides the spin susceptibility, XspmiT), temperature independent 
diamagnetic orbital contributions, Xdia; from the core electronic shells and the Van Vleck 
paramagnetism, Xvv? from the open shells of the Cu^"^ ions become important. Xdia can be 
estimated from the increments using Selwood's tabulated values as -12x10^^ cm^/mole and 
-26x10^^ cm'^/mole for Cu^"*" and Cl~, respectively, giving a total temperature independent 
diamagnetic susceptibility for CuCl2 of -64xl0~® cm^/mole.— An estimate of the Van Vleck 
contribution, xvv can be obtained from (cf. Ref [l^. 



Xtip 



A ' 



(2) 



where A'a is Avogadro's constant and A = 10|-Dg| is the energy separation of the d 
orbital states in a cubic crystal field, which amounts to approximately 25x10'^ cm~^. With 
Np^^\/kB ~ 0.375 cm^ K/mole one estimates a Van Vleck susceptibility of Cu^"*" in a cubic 
crystal field of the order of -|-43x 10"^ cm^/mole. In CuCl2 due to the Jahn- Teller distortion 
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the symmetry is reduced and the Van Vleck contribution depends on the actual matrix 
elements and the energy separation of the orbital states to the dx2-y2 state. For a detailed 
discussion see for example Ref. [l5|. Since these energy splittings are not known at present, 
we used an alternative method to determine the temperature independent contributions 
Xo'- In order to perform this analysis, we assume that the magnetic susceptibility of CuCl2 
can be described by a sum of the temperature dependent spin susceptibility Xspm(T) and a 
temperature independent contribution xo 

X{T) = Xo + Xspm(T) (3) 

where Xo is composed by the sum of the diamagnetic and the Van Vleck contribution 
according to 

Xo = Xdia + Xvv- (4) 

At sufficiently high temperatures, T >> 9, the spin susceptibility Xspm(T) can be ap- 
proximated by a Curie-Weiss law according to 

Na9'^^IS{S + 1) 

When fitting the high temperature susceptibility to the sum of equations (jlj) and ([5]) it 
is found that the fitting parameters, g, Xo ^-nd ^-re strongly correlated, and meaningful 
results cannot be obtained without fixing one of them. In order to reduce the number of 
parameters we therefore determined the g-factor, independently by an EPR experiment 
on a crystal (see below) to g^ = 2.050(1). The EPR experiment also revealed that our CuCl2 
crystals consist of mirror twins in the ab crystallographic plane with the mirror twins having 
the same b axis. By using gi, = 2.05 and fitting Eq. ([3]) to the (uncorrected) susceptibility 
data for T > 150 (see Figure [5]), we determine the following parameters: 

Xo = -18(1) X 10-^cm7mole, 
6 = -71(1)K. 
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FIG. 5: (Color online) Reciprocal magnetic susceptibility (uncorrected for diamagnetic and Van 
Vlexck contributions) of a crystal of CuCl2 determined in an external magnetic field of 1 T oriented 
parallel to the crystallographic 6-axis. The (red) solid line represents a Curie- Weiss law including 
a temperature independent term xo fitted to the data with parameters as given in the text. 

B. Magnetization Measurements 

Laboratory based magnetization experiments at 1.8 K with -Bcxt || b evidenced a spin- 
flop transition near 4 T. Measurements in a pulsed magnetic field at 1.4 K up to ~ 50 T 
confirmed this result. Evidence for further magnetic phase transitions at fields larger than ~ 
4 T could could not be found. Figure [6] displays the magnetizations measured in the pulsed 
field experiment. Saturation is not seen up to 50 T. The 1.4 K magnetization trace shows a 
slight upward deviation from nonlinearity while the 4.2 K traces within experimental error 
show a linear dependence of the magnetization over the full field range. 

Magnetic susceptibility experiments in various magnetic fields reveal that Tn depends only 
weakly on the field. The combination of the magnetization and the magnetic susceptibility 
experiments allows us to construct a magnetic phase diagram for Bcxt \\ b a.s displayed in 
Figure[7]with a bicritical point at 5.4(1) T and 24.2(1) K where the paramagnetic, spin flop 
and Neel phases coexist. 
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FIG. 6: (Color online) Magnetization of a crystal of CuCl2 measured at 1.4 K and 4.2 K in a pulsed 
magnetic field with the magnetic field aligned along b. The inset shows a magnification of the low 
field data. The solid and dashed lines represent data at 1.4 K and 4.2 K, respectively, measured 
with increasing (solid line) and with decreasing (dashed line) magnetic field. 

C. Electron Paramagnetic Resonance Measurements 

In order to determine the (^-factors along the principal axes of the (yf-tensor we carried 
out EPR measurements on small crystals of CuCl2 at room temperature at a microwave 
frequency of 9.46 GHz (X-band) (resonance field ~ 0.33 T). In a monoclinic system, one 
principal axis of the gf-tensor is fixed by symmetry and oriented parallel to the twofold axis, 
i.e. the crystallographic b axis for CuCl2. The remaining two principal axes of the ^f-tensor 
lie within the ac plane. To obtain the (^-factor along b, a crystal was oriented with the b-axis 
within the magnetic field plane. Rotation was carried out around an arbitrary axis perpen- 
dicular to b. A single resonance line is observed. Figure [8] shows the angular dependence 
of the (yf-factor obtained from fits of a field derivative of a Lorentzian to the resonance line. 
The angular dependence of the (7— factor shows the typical oscillating behavior. The minima 
occur when the field is oriented parallel to the b axis. The ^f-factor with field oriented parallel 
to the 6-axis amounts to 



The linewidths of the resonance lines exhibit angular dependence and they vary between 



gb = 2.050(1). 
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FIG. 7: (Color online) Magnetic phase diagram of CuCl2 derived from magnetization and heat 
capacity measurements with the magnetic field aligned along b. The inset displays the derivative 
dM/dH for selected fields H from which the phase boundaries have been extracted. 
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FIG. 8: (Color online) Angular dependence of the g'-factor of the EPR resonance line of Cu^"*" 
in a crystal of CuCl2 when the crystallographic 6-axis lies within the magnetic field plane. With 
b II i?ext the g-factor amounts to 2.050(1). 

20 mT and 40 niT due to small anisotropic exchange components.- Details will be discussed 
in a forthcoming paper.— 

If the magnetic field was oriented in the ac plane of the crystal (rotation around 6-axis), 
we observed an EPR spectrum that consists of two overlapping resonance lines, indicating 
that the crystals were twinned. The resonance fields and the linewidths of each line and the 
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FIG. 9: (Color online) Symbols: EPR spectra of CuCl2 with the external magnetic field oriented 
in the crystallographic ac-plane at (arbitrary) angles as given in the inset. The solid lines are fits 
to the measured spectra assuming two overlapping derivatives of Lorentzian resonance lines with 
separate resonance fields, linewidths and intensities, respectively. 

intensity were obtained by fitting to the spectra the sum of two Lorentzian absorption hnes 
(see Figure [9]). 

The (^-factors of the two hnes are shown in Figure [TOl One observes the typical oscillating 
angular dependence of the (^-factors on the angle between magnetic field the crystal orien- 
tation. Both lines show the same oscillating angular dependence of the (jf-factor, however 
with respect to each other shifted by ^57°. The extremal values are identical within exper- 
imental errors. The intensities for the two crystallites are about the same indicating that 
the investigated twinned crystallites have about the same volume. At angles where both 
twin-individua have the same resonance field and (^-factors, it is difficult to deconvolute the 
two resonance line unambiguously. This leads to the rather high error level in the fitted 
parameters especially at angles where the two gf-factor curves intersect. 

The principal ^f-factors, gi and g2 were determined by comparing the results with the 
general angular dependence of the ^f-factor for an axial system (Ref. Il6l ) 



where is the angle between the external field and the direction of the principal axis 1. 




(6) 
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FIG. 10: (Color online) Symbols: EPR spectra of CuCl2 with the external magnetic field oriented 
in the crystallographic ac-plane at (arbitrary) angles as given in the inset. The solid lines are 
fits to the measured spectra assuming two overlapping Lorentzian resonance lines with individual 
resonance fields, linewidths and intensities each. 

00 is a phase factor adjusted individually for each individual crystal. 
The solid lines in Figure [10] were calculated according to Eq. ([6]) using 



This valueis in good agreement with the ^f-factor obtained from the powder susceptibility 
measurement (see. Figure Hj). 

g2 and gb are almost identical indicating that the local symmetry for Cu^"'' in CuCl2 is 
very close to axial with gi = g\\ =2.255 and gt ^ 92 = 9± ^ 2.05. Here, we have identified 
the large gi with the (^-factor along the quasi-fourfold symmetry axis, perpendicular to the 
CuCU squares in the ribbons and gb and g2 with the ^f-factors for the magnetic field aligned 
within the planes of the ribbon. 



gi = 2.255(2) 



g2 = 2.055(2). 



With gb = 2.050(1), these values result in an average (^-factor of 




(7) 



g = 2.120(5). 
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This finding reflects the local symmetry, which is very close to tetragonal given the 
elongated CuCle octahedron containing each Cu^"*" ion. The four CI atoms in the equatorial 
plane are all at the same distance of 2.253 A and the apical CI atoms at a distance of 2.911 A 
(T = 42.9 K); the normal of the equatorial plane and the direction to the apical CI atoms 
encompass an angle of only 3.2°. 

The deviation of the (^-factors of the ground state of a 3(f conflguration in octahedral 
fleld with tetragonal splitting from the g-ia.ctoT of the free electron is determined by the 
spin-orbit coupling parameter A = -(/2S, where ( is the one-electron spin-orbit coupling 



parameter and S'=l/2. For Cu^"*" ( amounts to 829 cm ^ (Ref. Il6l ) and the deviations from 
the free electron ^f-factor are given by 

8A 

511-2 = ^ 

^0 

2A 

(8) 

where Aq and Ai are the energy separations between the d^2_y2 state and the d^y and the 
dyz,xz electronic states, respectively.— Using our experimental results, Aq and Ai amount 
to ~ 26x10^ cm~^ and ~ 30x10^ cm~^, respectively, consistent with the estimate we made 
from the magnitude of the temperature independent Van Vleck contribution. 

The principal axes of the two crystallites are shifted by 180° - /5 ~ 57° with respect to 
each other, where f3 is the monoclinic angle (cf. Table [T]). This flnding indicates that CuCl2 
crystals generally are twinned with the twins related to each other by a mirror symmetry 
operation with the mirror symmetry plane in the ac plane making the crystals to appear 
pseudo-hexagonal. This flnding becomes essential for the assignment of the propagation 
vector of the magnetic structure from elastic neutron diffraction on CuCl2 crystals. 

V. THERMAL PROPERTIES 

A. Heat Capacity Measurements 



The heat capacity of anhydrous CuCl2 has been determined before by Stout and Chisholm 
(Ref. ^ and by Billerey et al. (Ref. [l^. Figure [TT] shows our heat capacities Cp of a CuCl2 
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crystal and a powder sample. Our data agree very well with the data published by Stout 
et al. and those of Billerey and collaborators. The data of Stout and Chisholm are shown 
in Figure [TT] for comparison. The heat capacities of our powder sample and the crystal 
are very similar, and they exhibit A-shaped anomalies due to the onset of long-range AFM 
ordering. The anomalies are centered at ~23.8 K (crystal: 23.87 K; powder 23.82 K) in 
good agreement with the values reported by Stout et al. (23.91±0.1 K) and Billerey et 
al. (23.9±0.1 K). The maximum value in the A-anomaly is somewhat larger than the value 
reported by Stout et al. and Billerey et al. possibly indicating that our samples structurally 
are better ordered and contain less structural defects. The A-anomaly does not change its 
shape nor shifts the position when a magnetic field of up to 9 T is applied perpendicular 
to the crystal h axis. Sufficiently below the anomaly where the lattice contribution can be 
approximated by a Debye-T^ power law, the total heat capacities follow well the same power 
law Cp = Cph+Cmagnon oc indicating magnon contributions from a three-dimensional 
antiferromagnet . 




FIG. 11: (Color online) Temperature dependence of the specific heat of CuCl2. The main frame 
shows the heat capacity of a crystal, the lower inset shows the magnetic field dependence of the 
anomaly at Tn = 23.85(3) K. The magnetic field was directed perpendicular to the b axis. The 
upper inset displays the heat capacity near Tn of a polycrystalline sample. 

The entropy removed by long-range ordering signaled by the anomaly at 23.8 K amounts 
to 0.60(5) J/molK equivalent to 10% of _Rln2 as expected for a S=j spin system. The 
remaining 90% of the entropy apparently are removed by short-range ordering above T^. 
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The total heat capacity, Cp, of CuCl2 contains the spin (magnetic) and contributions 
from the phonon system (lattice contribution) according to 



Cp(T)=C^ag(T) + aatt(T) (9) 

In order to separate the magnetic contributions, Cmag, from the total heat capacity of 
CuCl2, one has to determine the lattice heat contributions, Ciatt- The lattice contributions 
can be estimated by fitting a power series expansion to the high temperature regime, where 
magnetic contributions are negligible, and extending the series to low temperatures. A sig- 
nificantly improved estimate of the lattice heat capacity is usually achieved by measuring the 
heat capacities of isostructural nonmagnetic compounds. The lattice reference systems ap- 
propriate for CuCl2 would be MgCl2 and CdCl2. Both compounds crystalhze with a layered 
crystal structure similar to that of CuCl2.— However, the MClg (M=Mg, Cd) octahedra of 
MgCl2 and CdCl2 do not exhibit a Jahn- Teller distortion, and consequently the two-valent 
anions are coordinated regularly by six CI atoms at the same distance of 2.64 A and 2.96 A 
for CdCl2 and MgCl2, respectively.—. 

Billerey et al. determined the heat capacities of MgCl2 in the temperature range 2 K 
< T < 100 K and constructed the lattice contribution to the heat capacity of CuCl2. We 
re-measured the heat capacities of a polycrystalline sample and of crystalline pieces of high 
purity (Alfa Aesar 99.99%) MgCl2 in an extended temperature range. We found that our 
data for MgCl2 obtained on two independent samples with two different experimental setups, 
especially, in the temperature regime relevant for the magnetic part below 100 K deviate 
markedly from Billerey's data. Our values being up to 10% smaller then Billerey's results. 

Additionally we determined the heat capacity of a small polycrystalline beads, typical 
diameter 2mm) of CdCl2 (Alfa Aesar 99.998%) The heat capacity of CdCl2 has been mea- 



sured before by Itskevich et a/.(Ref. |20| ) in a very limited temperature range. Takeda et 
al. (Ref. I21I ) also determined the heat capacity of CdCl2 in order to subtract the lattice 
contribution to the heat capacity of magnetic VCI2, but no data are available from the work 
of Takeda et al. Only the low-temperature cubic term oc is quoted. Figure [T^ shows the 
heat capacities of MgCl2 and CdCl2 as determined in this work in comparison with literature 
data. 

Figure [12] displays the heat capacities of CUCI2 and of MgCl2 and CdCl2 for comparison. 
As expected, the heat capacity of MgCl2 is significantly smaller, that of CdCl2 markedly 

18 




r(K) 



FIG. 12: (Color online) Temperature dependence of the specific heat of a single crystal of 
(black o) and of CdCl2 (blue □) compared with literature data as quoted in the inset (Refs. 
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20). 



greater than that of CuCl2. A detailed discussion how to construct an adequate lattice 
reference for CuCl2 and the magnetic contributions to the heat capacity of CuCl2 will be 
given in Section rVIIIi 



VI. NEUTRON DIFFRACTION 



A. Elastic Magnetic Scattering of a Powder Sample and a Twinned Crystal 

In order to search for magnetic Bragg reflections of the long-range ordered AFM state we 
examined in detail the low-angular parts of the powder diffraction patterns. A comparison of 
two diffraction patterns collected in D20's highest intensity mode at a wavelength of A=2.4 
A at 2 K and 30K reveals additional magnetic scattering up to 20 ~ 40° (Figure [T3](a)). 
Most pronounced is a reflection ~ 23° at the low angle shoulder of the nuclear 01 reflection. 
Additional, though significantly less intense, magnetic scattering intensity is visible in the 2 
K pattern at ~34°. Both features are not present any more in the pattern collected at 30K 
and thus can be clearly identified as magnetic Bragg reflections. 

Towards low scattering angles, the 30K pattern lies clearly above the 2 K pattern indi- 
cating significant paramagnetic scattering in the 30 K pattern. In the difference pattern the 
magnetic Bragg reflection near 19° is slightly deformed owing to a small shift of the 01 
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nuclear reflection arising from the contraction of the lattice between 30 K and 2 K. 




FIG. 13: (Color online) Neighborhood of the 001 nuclear Bragg reflection of CuCl2 (A=2.4 A ) at 
2 K and 30K. The additional magnetic Bragg reflection near 20 w 23.2° is clearly visible. It can 
be indexed as -1,0 225 0.5. 

In order to obtain the temperature dependence of this first magnetic Bragg refiection, we 
collected a series of patterns between 2 K and 45 K in the D20's higher resolution mode using 
a wavelength of A= 1.888 AThe difference between the pattern collected at the respective 
temperatures and the pattern at 28.9 K was fitted to a Gaussian the intensity of which is 
displayed versus temperature in Figure [141 

The temperature dependence of the integrated intensities can be well fitted by a critical 
power law according to Eq. (fTOj) with the critical exponent (3 

I{T)=h{l-T/T^fP. (10) 

A fit of the integrated intensities converges to a Neel temperature Tn=24.7(1) K consistent 
with the magnetic susceptibility and the heat capacity and a critical exponent /3=0.16(1) 
(see Figure [TW a)). 
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FIG. 14: (Color online) (a) Temperature dependence of the integrated intensity of the magnetic 
Bragg reflection observed at 26 ^ 18.2° (A=1.888 A ) (o). The sohd line is a fit of the experimental 
data with a power law with a critical temperature Tn=24.7(1)K. (b) Temperature dependence of 
the integrated intensity of the magnetic Bragg reflection collected on a crystal of CuCl2 (DIO four- 
circle diffractometer) A=2.36 A ) (o). The solid line is a fit of the experimental data with a power 
law with a critical temperature Tn=23.9(1)K and a critical exponent /3=0.23(2). (c) log-log plot 
of the integrated intensity displayed in (b) vs. the reduced temperature t=l-T-^/T. The solid line 



Intensities of a set of magnetic Bragg reflections were collected on a twinned crystal. 
Figure [T^b.c) display the intensity of a strong magnetic reflection collected on a twinned 
crystal on ILL's four-circle diffractometer DIO.-^*^^ Initially, this reflection was indexed as 
0.776 0.5. In view of the magnetic structure refinement of the powder data (see below) we 
must conclude that this reflection originates from the related twin individuum in the crystal. 
Using the inverse of the twin matrix eq. (1) the indexes transform to 0-0.776-0.5 i.e. 110 
- r and 1 1 1 - r , where r is the magnetic propagation vector refined from the powder data 
below. A fit of the integrated intensities of this reflection (solid line in Figure [TWb.c)) gives 
Tn=23.9(1) and /9=0.23(1) consistent with the powder diffraction results. 

Although the investigated temperature regime is limited, and corrections for critical scat- 
tering near Tn have not been applied, the refined critical exponents (3 obtained from the 
powder and crystal are significantly lower than those characterizing the well known three- 
dimensional universality classes (Ising, XY and Heisenberg) for which (5 ranges between 
0.325 and 0.365.— The refined critical exponent (3 ~ 0.2 for CuCl2, however, comes close to 
the predictions of Kawamura for the universality class of of the SO (3) or the Z2 x Si sym- 
metry for systems with chiral order. For such systems (3 ranges between 0.22 and 0.25.—'^'^ 
Figure \\Mc) displays the intensity of a magnetic reflection collected on a twinned crystal on 
ill's four-circle diffractometer DIO together with a power law with a critical exponent (3 = 
0.24(2) 

The difference between the 2 K and the 30 K powder diffraction patterns (lower panel of 
Figure [T3i) reveals an additional weak triple of magnetic Bragg reflections between 32° and 
37° with a peak shape very similar to those of neighboring nuclear reflections. Additionally, 
the difference patterns allows one to identify magnetic scattering around a scattering angle 
of -70°. 

All magnetic Bragg reflections can be indexed with a propagation vector (1, 0.225,0.5) 
indicating a doubling of the magnetic cell along the c axis and an incommensurate AFM 
ordering along h (see Figure fT5l). 

The magnetic structure of CUCI2 was refined using the D20 2 K and the difference powder 
diffraction data sets using the program Fullprof. In order to obtain reliable parameters 
for the nuclear structure (nuclear scale factor, profile parameters), we first refined the 2 
K pattern assuming nuclear scattering only. Since magnetic scattering is very small this 
refinement converges well and provides the lattice and profile parameters and the scale 
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FIG. 15: (Color online) Neutron diffraction pattern (difference of 2 K and 30K patterns) (red o) 
(A=2.4 A ) of a powder sample of CuCl2 collected at the high intensity diffractometer D20 plotted 
in comparison with the pattern obtained from a profile refinement assuming the magnetic structure 
described in detail in the text (solid black line). The difference is shown in the lower part by the 
(blue) solid line. The positions of the the magnetic Bragg reflections used to simulate the patterns 
are indicated by the small (green) vertical bars in the lower part of the figure. 

factor. Subsequently a difference pattern obtained by subtracting the 30 K pattern from 
the 2 K pattern was used to refine the magnetic structure on the basis of the propagation 
vector and the nuclear scale factor, which was kept fixed in the refinements. The background 
of the difference pattern was chosen manually and subtracted. The refinement based on 18 
magnetic reflections was successful assuming the nuclear space group CI and a helix with the 
moments confined to the be plane. Several tests were carried with the moment directions 
pointing in a more general direction out of the be plane but did not lead to significantly 
improved fits. 

Assuming this model for the magnetic structure the ordered moment at 2 K converged 

to 



/i(2K) = 0.50(l)/iB. 
The propagation vector r at 2 K was refined to 

f(2K) = (1,0.2257(6), 0.5) 
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implying an angle of 81.2(3)° between neighboring moments along the chain. Figure [16] 
displays the magnetic structure of CuCl2 at 2 K. 




FIG. 16: (Color online) Magnetic structure of CuCl2 at 2 K. A helix propagating along b with the 
moments confined to the be plane. Neighboring moments along b enclose an angle of ~ 81°. 



VII. DENSITY FUNCTIONAL CHARACTERIZATION OF THE SPIN EX- 
CHANGE INTERACTIONS AND THE SPIN-SPIRAL STATE 

To discuss the magnetic structure of CuCl2, we consider the intra-chain and the inter- 
chain spin exchange interactions defined in Figure [T71 Ji and J2 are the NN and NNN 
intra-chain spin interactions, respectively. J3 and J4 are the inter-chain spin exchange inter- 
actions along the c-direction, and J5 is that along the a-direction. The structural parameters 
associated with the spin exchange paths Ji - J5 are summarized in Table [Tll 

To estimate the spin exchange parameters Ji - J5, we first determine the total energies 
of the six ordered spin states of CUCI2 presented in Figure [18] on the basis of first principles 
density functional theory (DFT) electronic band structure calculations. Our DFT calcula- 



tions employed the Vienna ab initio simulation package (Ref. l26lJ27lJ28l ) with the projected 



augmented-wave method, the generalized gradient approximation (GGA) for the exchange 



and correlation functional, (Ref. l29l ) the plane-wave cut-off energy of 330 eV, and the sam- 
pling of the irreducible Brillouin zone with 36 fc-points. To take into consideration of the 
strong electron correlation associated the Cu 3d state, we performed GGA plus onsite repul- 
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FIG. 17: (Color online) Definition of the five exchange paths Ji - J5 of CuCl2, where the numbers 
1-5 represent Ji, ... J5, respectively. 



path 



Cu...Cu (A 



C1...C1 (A 



CU-C1...C1 (°) 



Ji 



3.299 



J2 



J3 



6.599 



6.824 



3.299 (x2) 
3.729 (x2) 



136.^ 



133.1 



9.493 



4.979 



174.1 



J5 



7.608 



3.808 



131.7 



TABLE II: Geometrical parameters associated with the spin exchange paths Ji - J5 of CuCl2. 



sion (GGA+U) calculations (Ref. l30l ) with f/ = 3, 5 and 7 eV on Cu. The relative energies 
of the six ordered spin states obtained from our GGA+U calculations are summarized in 
Table [ml 

The energies of the six ordered spin states can also be written in terms of the spin 
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state 



U = 3eV 



U = 5eV 



U = 7 eY 



FM 



AFMl 



-11.23 



-8.62 



-6.52 



AFM2 



21.68 



20.06 



17.43 



AFM3 



13.77 



14.09 



12.95 



AFM4 



-31.01 



-22.12 



-15.93 



AFM5 



-39.72 



-28.85 



-21.05 



TABLE III: Relative energies (in meV) of the five ordered spin states of CuCl2 obatined from the 
GGA-I-U calculations. 



FM 





AFl 



AF2 




AF4 



FIG. 18: (Color online) Six ordered spin arrangements of CuCl2 employed to extract the spin 
exchange parameters Ji - J5. The up- and down-spin Cu sites are represented by the open and 
filled circles, respectively. 
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Hamiltonian, 



H = -Y,J^J^A^ (11) 
(1) where (= Ji - J5) is the spin exchange parameter for the spin exchange interaction 
between the spin sites i and j, while Si and Sj are the spin angular momentum operators 



at the spin sites i and j, respectively. By applying the energy expressions obtained 
dimers with N unpaired spins per spin site (in the present case, N = 1), (Ref. l31 



or spin 



32) the 



total spin exchange energies per formula unit of the six ordered spin states (see Figure [T5]) 
are written as 





(-2Ji 


-2J2 


-2J3 


-4J4 


-4J5) 


Ar2 
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Eafi = 


(-2Ji 


-2J2 


+ 2J3 


+ 4J4 


+ 4J5) 
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EaF2 = 


(+2Ji 


-2J2 


-2J3 


-4J4 


+ 4J5) 
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EaF3 = 


(+2Ji 


-2J2 


+ 2J3 


+ 4J4 


-4J5) 




Eaf4 = 


(+2Ji 


-2J3 


+ 4J4 








EaF5 = 


(+2Ji 


+ 2J3 


-4J4 









(12) 

By mapping the relative energies of the six ordered spin states given in terms of the 
spin exchange parameters to the corresponding relative energies obtained from the GGA+U 
calculations, we obtain the values of the spin exchange parameters Ji - J5 summarized in 
Table lYl 

For all values of U employed, the two strongest spin exchange interactions are the intra- 
chain interactions Ji and J2. As found for the magnetic oxides with CUO2 ribbon chains, 

(Ref. yO) Ji is ferromagnetic while J2 is antiferromagnetic with J2 larger in magnitude 

than Ji. As a consequence, the intra-chain spin exchange interactions are geometrically 
frustrated. Except for J3, the inter-chain spin exchange interactions are negligible; J3 is 
antiferromagnetic and is weaker than J2 by a factor of approximately four. Because of the 
intra-chain spin frustration, the magnetic ground state of each CuCl2 chain is expected to be 
a spin-spiral state with a certain incommensurate repeat vector qy along the chain direction, 
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state 



U = 3eV 



U = 5eV 



U = 7 eY 



Jl 23.3 21.4 18.4 

Ja -41.4 -31.9 -24.5 

J3 -9.1 -7.0 -5.3 

J4 -0.2 -0.1 -0.1 

J5 -0.8 -0.7 -0.5 

TABLE IV: Values of the spin exchange parameters Ji - J5 (in meV) in CuCl2 obtained from the 
GGA-I-U calculations. 

i.e., the 6-direction. The inter-chain spin arrangement is expected to be antiferromagnetic 
along the c-direction due to the spin exchange J3. Our calculations of the energy E (0, Qy, 
1/2) as a function of the modulation wave vector q = (0, qy, 1/2) shows a minimum at qy ~ 
0.22. 

The magnetic state with the spin-spiral arrangement removes inversion symmetry and 



hence should induce ferroelectric polarization. 



b confirm this prediction, we carried out 



GGA-I-U plus spin-orbit coupling (SOC) (Ref. l35l ) calculations for CUCI2 that has the com- 
mensurate spin-spiral arrangement with q = (0, 0.25, 0) to reduce the computational task. 
Indeed, our calculations of the ferroelectric polarization using the Berry phase method (Ref. 
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371 ) for the spin-spiral state of CnCh determined by the GGA-fU-(-SOC calculations lead 
to nonzero ferroelectric polarizations; with the spin-spiral plane parallel (||) and perpendic- 
ular (± ) to the CUCI4 ribbon plane, the x-, y- and z- components of the polarizations are 
calculated to be P\\ = (-44.5, 0, 71.6) and = ( 5.8, 0, 28.3) in units of /iC/m^. 
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VIII. EXPERIMENT VERSUS THEORETICAL RESULTS 



In the following we will compare the implications of the theoretical results for the magnetic 
susceptibility and the thermal properties. The calculations indicate that the intra-chain 
exchange is dominant and that the magnetic properties of CuCl2 have to be described as 
those of a 5*= 1/2 frustrated Heisenberg chain with FM nearest and AFM NN interaction. 
According to the calculations, the ratio J1/J2 is negative and is expected to be of the 
order of ~ -0.5 to -0.75. Additionally, there is small but noticeable inter-chain interaction 
being dominant along the c axis that leads to cooperative AFM ordering and the AFM 
arrangement along the c direction with the doubling of the magnetic unit cell along c. 
Magnetic and thermal properties of S=l/2 zig-zag Heisenberg chains have been calculated 
in great detail by exact diagonalization by Heidrich-Meisner and coworkers, and numerical 
tables are available for a wide range of the ratio a = Ji/ J2-—'^ These works comprise also 
the magnetic susceptibility of a S=l/2 Heisenberg chain with NN AFM exchange only, i.e., 
a=0 as calculated before by Kliimper et al. and parametrized by a Fade approximant by 
Johnston et al..—^ 



A. Magnetic Susceptibility 

Correcting the magnetic susceptibility results with the sum of the temperature indepen- 
dent diamagnetic and Van Vleck corrections to the susceptibility (xo) gives the pure spin 
susceptibility Xspin(7') of CuCl2. Using the (7-factor along b as determined by the EPR mea- 
surements, (7^=2.050(1), enables us to compare the magnetic susceptibility of a CuCl2 crystal 
measured with magnetic field along the 6-axis with model calculations. The only remaining 
adjustable parameters are the intra-chain exchange constants Ji and J2- The effect of the 
finite inter-chain exchange Jintcr on the paramagnetic susceptibilities calculated for the zig- 
zag chains has been taken into account by employing a mean field approach, which gives as 



inter-chain interaction corrected susceptibility as(Ref. |40| ) 



cor /rp\ _ AspinV-t ) 

Aspinl^ ; T T/M. ,,2\ frr\- > 



Xspin(^) 
l-(^J/iVAS4)Xspin(T)' 

The leading intra-chain exchange J±_ is provided by J3 (see Table UTTl and Figure [Hj) and 
amounts to about 20% of the next-nearest exchange J2 along the chain. Including intra- 
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FIG. 19: (Color online) Temperature dependence of the magnetic susceptibility of a crystal of 
CuCl2 measured with magnetic field aligned along the crystal 6-axis corrected for the sum of dia- 
magnetic and Van Vleck contributions (xo) (°) compared with the results of exact diagonalization 
calculations of 5=1/2 Heisenberg chains with NN and NNN exchange Ji and J2, respectively defin- 
ing the ratios a=Ji/ J2 as given in the inset. The (colored) solid lines were calculated according 



to Eq. m with 

Xspin(^) taken from Ref. 
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39l for J2= -117 K and as as given in the inset. Prom 



top to bottom Q amounts to -0.75, -1-0.25, respectively. 

chain exchange according to Eq. [13] renormalizes the calculated susceptibilities by about 8% 
near the maximum at ~73 K. At room temperature the decrease amounts to ~4%. 

Figure [19] shows the comparison of the spin susceptibility of CuCl2 as determined in this 
work with the inter-chain exchange corrected susceptibility of a Heisenberg chain with NN 
and NNN exchange interaction Ji and J2, respectively.— 1^ 

According to this comparison the magnetic susceptibility of CUCI2 is best reproduced by 
an inter-chain exchange corrected susceptibility of a Heisenberg chain with NNN exchange 
interaction J2=-117 K and a ratio a ~ -0.25, i.e., ferromagnetic nearest neighbor exchange 
as predicted by the DFT calculations. However, the ratio a is markedly smaller than pre- 
dicted by the DFT calculations (see Table [In|) . Including inter-chain exchange interaction 
is essential to reproduce the measured susceptibility correctly. 
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B. Magnetic Heat Capacity 



The magnetic contributions to the heat capacity of CuCl2 has been evaluated before 
by Stout and Chisholm and subsequently by Billerey and coworker s.-^^^ Our measurements 
confirm (see Figure [TTT a)) their peak-anomaly at ~23 K indicating long-range AFM ordering. 
Our data agree very well with those of the previous works, and the anomaly is somewhat 
sharper in our experiment possibly due to the fact that polycrystalline samples have been 
used beforc-ii^ Stout and Chishom and Billerey et al. used the heat capacities of MnCl2 and 
MgCl2, respectively, in order to construct an estimate of the lattice heat capacity of CuCl2. 
We redetermined the heat capacities of MgCl2 and found marked deviations of our data 
from those of Billerey. We considered also CdCl2 as another possible system to estimate the 
phonon contributions to the heat capacity of CuCl2 and measured the heat capacity of CdCl2 
over a temperature range much wider than covered by the experiments of Itskevich et al.{see 
Figure [TWa)). MgCl2 and CdCl2 are diamagnetic and have a crystal structure containing 
MgCl2 and CdCl2 layers, which closely resemble those in the crystal structure of CuCl2. 
However, the MClg (M=Mg, Cd) octahedra are not axially elongated as they are in the CuOg 
octahedra. In addition, owing to the different atomic masses of Cu and Mg, the phonon 
spectrum of CuCl2 and MgCl2 and Cd2 may be expected to be different, especially at low 
energies. Although MnCl2 is magnetic, its AFM ordering occurs at very low temperature, 
1.81 K and 1.96 K, and magnetic contributions to its heat capacity are negligible above 
~10K.— 1^ Above lOK, MnCl2 may therefore provide a reasonable approximation to the 
lattice heat capacity of CuCl2. In order to compensate for the difference in the cation masses, 
we employed the method of corresponding states and assume that the phonon spectra of 
MgCl2, CdCl2 and MnCl2 are very similar to that of CuCl2 and the different phonon energies 
can be compensated for by modifying the temperatures.— Compensation of the different 
cation mass was done by stretching or compressing the temperature scale by a constant 
factor. This factor was chosen such that at sufficiently high temperatures where magnetic 
contributions to the heat capacity vanish the heat capacities of the lattice references match 
with those of CuCl2 (see Figure fTT]) . We found factors of 0.88, 0.93 and 1.18 to be appropriate 
for MgCl2, MnCl2 and CdCl2, respectively. 

A comparison of the heat capacities of CuCl2 and MCI2 (M=Mg, Mn, Cd) reveals sig- 
nificant short-range ordering magnetic contributions to the heat capacity of CuCl2 up to ~ 
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100 K, as already concluded by Stout et al. and Billerey^ii^. 

In order to construct a lattice heat capacity reference for CuCl2 from the heat capacities 
of MgCl2, MnCl2 and CdCl2, we tried various combinations of the heat capacities of MgCl2, 
MnCl2 and CdCl2 with temperatures scaled by the factors given above and subtracted these 
from the total heat capacity of CuCl2. The resulting magnetic heat capacity of CuCl2, 
Cmag/r was subsequently integrated and the total magnetic entropy was compared with 
Mn2, the entropy expected for a 5=1/2 magnetic system as a crosscheck. 

Due to the very large deviations of the heat capacity of MgCl2 from that of CuCl2, 
especially in the temperature regime below ~ 50K where magnetic contributions are essential 
we discarded MgCl2 and achieved an adequate lattice heat capacity reference by averaging 
for each temperature the heat capacities of MnCl2 and CdCl2. 

The magnetic heat capacity of CuCl2 obtained after subtracting the phonon contribu- 
tions is shown in Figure [20r b) and (c). In addition to the A-type anomaly at ~23 K the 
magnetic contributions to the heat capacity are characterized by a broad anomaly due to 
short-range AFM ordering with maximum at 35 K (in Cmag/T" and a maximum value of 
0.11 J/molK^). Our data are close to those of Billerey et al. but deviate markedly from the 
data published by Stout et al., who found a significantly higher value for the short-range 
ordering maximum. The total entropy of the magnetic heat capacity amounts to 0.66(2)i?, 
very close to B\n2 expected for the entropy of a spin system with S'=l/2. Figure [20r c) shows 
a comparison of the magnetic heat capacity with various model calculations, e.g., the heat 
capacity of an AFM Heisenberg chain with NN only and with NN and NNN interactions 
with exchange parameters Ji and the ratio J2I J\ similar to those discussed for the mag- 
netic susceptibility.—'^ Both chain models deviate significantly from the experimental data. 
Surprisingly good agreement of the broad magnetic heat capacity anomaly is obtained for a 
simple model assuming a spin 5*= 1/2 dimer according to 



i/=-Jdim^l^2, (14) 



with the heat capacity discussed in detail in Ref. |42|for an exchange constant Jdim ~ -117 
K. The data are also well reproduced by an anisotropic three-dimensional Heisenberg model 
on a cubic lattice with in-plane exchange Jy and interplane coupling J^.— With Jy ^ -70K 
and ~ J|| maximum value and position of the short-range ordering anomaly are rather 
well reproduced. The onset of three-dimensional ordering showing up as a well separated 
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FIG. 20: (Color online) (a) Temperature dependence of the specific heat of a single crystal of 
CuCl2 (black o) and of a polycrystalline sample measured by Billerey et al. (red □).— The heat 



44 ) are also given. 



capacities of MgCl2 and CdCl2 as measured in this work and of MnCl2 (Ref. 
Their temperature axes have been rescaled by the factor 0.88, 0.92 and 1.18 for MgCl2, MnCl2 
and CdCl2, respectively. The (black) solid line is the lattice reference (phonon contribution to the 
heat capacity) constructed from the heat capacities of MnCl2 and CdCl2 as described in detail in 
the text, (b) Magnetic contribution to the heat capacity of CuCl2 obtained as the difference of the 
CuCl2 data and the lattice reference displayed in (a). Published data by Stout et al. and Billerey 
et al. are displayed for comparison, (c) Magnetic heat capacity of CuCl2 (black o) compared with 
various model calculations as labeled in the inset and discussed in detail in the text. 
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peak well below the short-range ordering maximum appears close to the cooperative ordering 
temperature observed for CuCl2. 

IX. DISCUSSION 

Complex compounds containing Cu^"^ ions with low- dimensional magnetic behavior are 
legion and have been widely investigated.— However, the detailed magnetic behavior of 
anhydrous Copper (Il)-chloride, CuCl2, itself, however, has largely remained a mystery until 
now. From the early susceptibility studies by deHaas and Gorter and Starr et al. it had 
become evident that CuCl2 shows unusual AFM properties with a broad maximum in the 
magnetic susceptibility.-i^ The specific heat of CuCl2 measured by Stout and Chisholm 
proved long-range AFM ordering at around 23 K and the excess heat capacity extending 
well above was analyzed in terms of linear chain behavior. Because of a lack of better 
theoretical foundation this analysis was done using an Ising-type approach.- A neutron 
diffraction study in order to determine the long-range ordered groundstate has not been 
carried out so far. 

Our study presents a detailed confirmation of the early susceptibility study carried out on 
polycrystalline samples and twinned single-crystals of anhydrous CuCl2. Susceptibility and 
heat capacity confirm the one-dimensional AFM properties of CuCl2 and the Neel temper- 
ature 23 K as already found by Stout and Chisholm.- This rather high Neel temperature as 
compared to the maximum temperature of the short-range ordering susceptibility maximum 
occurring at about 70K indicates substantial inter-chain exchange interaction and underlines 
the conclusion by deJongh and Miedema that CuCl2 is a rather poor example of a chain 
structure.— 

CuCl2 consists of ribbons of edge-sharing CuC^ square planes with Cu^+ ions. In such 
systems a spin helicoidal groundstate can set in along the chain direction at low temperature 
due to the competition between NN ferromagnetic and NNN antiferromagnetic interactions 
along the Cu^^ chains. The first clear example proven to show such a behavior was LiCuV04 
that crystallizes with an inverse spinel structure and contains Cu02 ribbon chains. Helicoidal 
magnetic ordering with the helix propagating along the chain direction was found by a 
single-crystal neutron diffraction study.— The frustration scenario with ferromagnetic NN 
and larger AFM NNN exchange was proven by an inelastic neutron scattering study.— 
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Special interest in the magnetic properties of LiCuV04 arose lately from the observation 
of multiferroic behavior found by Naito et al. and Schrettle et al..—'^^ LiCuV04 shows 
a dielectric polarization below the Neel temperature and the polarization can favorably be 
switched by applying a magnetic field. The origin of this dielectric polarization is due to 
spin-orbit coupling on the Cu sites, however, the asymmetric charge density distribution 
necessary for the electric polarization occurs around the ions.— 

LiCu202 is another Cu^^ system with Cu02 ribbon chains exhibiting helicoidal spin 
order.— i^i^i^i^'^'^ In LiCu202, in contrast to LiCuV04, the magnitude of the FM NN 
(-11±3 meV) exchange is larger than the AFM NNN exchange (7±1 meV). Ferroelectricity 
in LiCu202 has been shown by Park and coworkers.— 

Another related system isostructural with LiCu202, which lately has attracted special 
attention, is NaCu202. A helicoidal magnetic ground state has been found by neutron pow- 
der diffraction and by NMR experiments.— >^>^>^ A review of recent systems with helicoidal 
ground states has been compiled by Drechsler and collaborators.— 

All helicoidal systems investigated so far contain oxygen ions to provide superexchange 
between the Cu^^ moments. Clear proof of the helicoidal ground state by magnetic neutron 
scattering is available only for a few of them. Evidence for systems containing anions other 
than O^"*" has not been reported yet. 

Anhydrous CuCl2 is therefore the first example for a halide system for which a helicoidal 
ground state has been proven unequivocally. Our electronic structure calculations give 
clear theoretical evidence for a scenario of competing FM NN and AFM NNN exchange 
interactions. Our experimental results for the bulk magnetic properties and the comparison 
with their predictions for a 5*= 1/2 AFM Heisenberg chain with competing FM NN and 
AFM NNN exchange support the theoretical predictions. Clear proof for the helicoidal 
AFM ground state comes from the neutron diffraction work on powder and single crystals. 
We find a magnetic structure indicating a S=l/2 helix in the CuCl2 ribbon chains. The 
component of the magnetic propagation vector along the &-axis amounts to 0.2257(6) and 
is in best agreement with the result of the electronic structure calculations. The spin-spiral 
as refined from diffraction data to is found to lie in the be plane, consistent with the spin- 
flop observed in the magnetization data when the magnetic field is oriented along the b 
axis. The refined magnetic structure indicates a magnetic unit cell which is close to being 
quadrupled as compared to the nuclear cell with neighboring moments enclosing an angle 
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of ~81° which leads to a significant reduction of NN exchange energy. With respect to the 
hehx this scenario is analogous to that in LiCuV4 (propagation vector (0,0.468,0) (Ref. l66l ) 
and NaCu202 (propagation vector (0.5,0.227,0).— 

The magnetic susceptibility data indicate the importance of inter-chain interaction, which 
according to the electronic structure calculation amounts to ~20% - 30% of the intra-chain 
interaction. Interchain interaction decreases the temperature for the short-range order max- 
imum and provides a better agreement of experimental data and theoretical prediction (Fig- 
ure [19]). 

Particularly surprising and striking are the deviations of the short-range ordering contri- 
butions to the magnetic heat capacity for CuCl2 and the heat capacity expected for a 5*= 1/2 
quantum chain with FM NN and AFM NNN exchange interactions (see Fig. I^UT c)). The 
short-range order contributions match significantly better with a 5*= 1/2 dimer singlet-triplet 
excitation scenario with AFM exchange ~ -117 K corresponding to the exchange found for 
the NNN interactions along the ribbons. We believe that this coincidence is due to the 
competing intra-chain competing interactions and the additional AFM inter-chain coupling. 
The magnetic configuration energy f/mag measures the short-range pair correlation function 
according to (cf. e.g. Ref. Il2l ) 



f/mag OC J( 5*0 5*1), (15) 

where the subscript 1 denotes a spin adjacent to 0. Thus f/mag and the magnetic heat 
capacity dUmag/dT in a frustrated chain system with coupling to neigboring chains may 
primarily sense magnetic correlations in a short-ranged cluster in which the configuration 
energy is determined by the leading AFM exchange to NNN spins. 

In summary, anhydrous CUCI2 shows ID AFM behavior and long range AFM ordering 
below a Neel temperature of 23.9 K, below which CUCI2 adopts an incommensurate magnetic 
structure (1,0.2257,0.5) with a spin-spiral propagating along b and the moments confined in 
the be crystallographic plane. The spin-spiral results from competing FM NN and AFM NNN 
spin exchange interactions along b. Anhydrous CUCI2 is the first halide quantum system 
containing CUCI2 ribbon chains for which a helicoidal magnetic ground state is realized, and 
is expected to be ferroelectric below 23.9 K. 
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